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We present the first ever calculation of cosmic microwave background (CMB) anisotropy power 
spectra from semilocal cosmic strings, obtained via simulations of a classical field theory. Semilo- 
cal strings are a type of non-topological defect arising in some models of inflation motivated by 
fundamental physics, and are thought to relax the constraints on the symmetry breaking scale as 
compared to models with (topological) cosmic strings. We derive constraints on the model param- 
eters, including the string tension parameter fi, from fits to cosmological data, and find that in this 
regard BPS semilocal strings resemble global textures more than topological strings. The observed 
microwave anisotropy at I = 10 is reproduced if Gfj, — 5.3 x 10 -6 (G is Newton's constant). However 
as with other defects the spectral shape does not match observations, and in models with inflationary 
perturbations plus semilocal strings the 95% confidence level upper bound is Gfi < 2.0 x 10 -6 when 
CMB data, Hubble Key Project and Big Bang Nucleosynthesis data are used (c.f. Gfi < 0.9 x 10 -6 
for cosmic strings) . We additionally carry out a Bayesian model comparison of several models with 
and without defects, showing models with defects are neither conclusively favoured nor disfavoured 
at present. 
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I. INTRODUCTION 



Recent observational data are tightly constraining cos- 
mological models, in many cases rendering them non- 
viable. One of the constraints comes from the fact that 
many potentially successful high-energy physics moti- 
vated models of inflation predict a cosmic string network 
[l|, l2| in the post-inflationary era: this occurs in super- 
symmetric hybrid inflation [3|, in grand unified theories 
4] and in some string theory inflation models [j| [|| 0, @] • 
These strings and other type of defects contribute to 
the cosmic microwave background (CMB) anisotropics 
in addition to the primordial inflationary fluctuations, 
allowing CMB experiments toprovide upper limits on 
the contribution from defects [3, [lCj, [n], [l2|, [l3| • In the 
near future, observations will either reveal the presence of 
these defects or will significantly tighten the constraints 
upon the models from which they are predicted (see for 
example Refs. [3 El El El El ) ■ 

However, a simple and elegant method for evading the 
constraints coming from cosmic strings in these mod- 
els was proposed in Refs. [19|, |2(|. A duplication of 
the matter fields responsible for the string formation, 
which arises quite naturally in several models [20J, [2l|, |22| , 
changes the nature of the strings from topological and 
therefore long-lived strings, to non-topological semilocal 
strings [23J ■ All the desired inflationary properties of the 
model persist, and semilocal strings are believed to be 
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cosmologically less harmful than their topological coun- 
terparts, in that their energy density may be lower for a 
given symmetry breaking scale. In this paper we confirm 
these expectations by performing the first CMB calcula- 
tions for semilocal strings and determine the constraints 
upon semilocal strings imposed by current data. 

Calculating the additional contribution that cosmic 
strings would make to the CMB temperature and po- 
larization power spectra is not an easy task. They are 
highly non-linear entities, and the range of scales in- 
volved in the problem is enormous. However, as the 
width of the string is much smaller than the string length, 
one possibility is to neglect the finite width and perform 
Nambu-Goto type simulations [H, H, [H, [27J]. Unfor- 
tunately, these simulations ignore the string decay into 
gravitational waves or particles and the corresponding 
back-reaction onto the network. Their representation of 
small-scale loop production in the network is also under 
debate [28|, [23, [30, [Ml, (23 • Yet another level of simplifi- 
cation often used in CMB calculations for local strings is 
to simulate merely a stochastic ensemble of unconnected 
string segments 0, El EI El, H Hi] . Lacking any dy- 
namical equations, in this case the sub-horizon decay of 
strings must be taken into account by the random re- 
moval of segments, at a rate chosen to match that seen 
in network simulations, but qualitatively accurate results 
are obtained with much less computational cost. In the 
global defect case, the energy is not localized into the 
string cores and their resolution is not essential. This en- 
ables the calculation of CMB contributions from global 
defects using field theory simulations [35|, (36| . 

None of these three approaches is useful for the study 
of semilocal strings, since they are not well-modelled by 
Nambu-Goto strings while, unlike the global case, the 
string cores are potentially important. Fortunately, in 
a recent breakthrough we performed the first CMB cal- 



culations for local cosmic strings to derive from field- 
theoretical simulations [33]. We achieved this by the use 
of a parametrized approximation to reduce the range of 
scales required in the simulations, but which had no re- 
solvable effect on the CMB predictions. That approach 
has opened the door for the possibility of simulating de- 
fects with no known approximating models. 

We have therefore employed the machinery developed 
for the local string case to the simplest model containing 
scmilocal strings, with the parameters chosen to give the 
Bogomol'nyi-Prasad-Sommerfield (BPS) limit (see Sec- 
tion [TTJ). The BPS limit arises naturally in D-term SUSY 
inflation and D3/D7 brane inflation. We leave the study 
of non-BPS cases to a future article but it is the expec- 
tation that in the two opposite extremes between which 
the BPS case lies, the semilocal result will tend to that 
of global textures or to that of local strings [38J, 123, [4fJ . 
We perform a detailed comparison of our results with 
these two cases using our earlier results in Ref. [33] for 
local strings. For the case of global textures, we perform 
CMB calculations using the linear a model. 1 We leave 
a discussion of our texture calculations to an Appendix. 

We have also performed a multi-parameter compari- 
son with current data, yielding results similar to those of 
Ref. [l3| where it was seen that local cosmic strings are 
mildly favoured by current CMB data. We find that the 
preference for semilocal strings is slightly greater and, as 
in Ref. [l3j, we explore this preference using Bayesian 
evidence calculations, as well as considering the implica- 
tions of non-CMB data. 

This paper is organized as follows: in the next Section 
we discuss the semilocal string model, before giving a 
brief overview of the simulation and CMB calculation 
method in Section IIIII We will then present our CMB 
results and consistency checks in Section IIV[ including 
the comparison with local strings and global textures. 
Section [V] then presents the comparison with data. 



II. SEMILOCAL STRING MODEL 

The simplest model giving rise to semilocal strings is 
given by the Lagrangian density 23] 
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The fields 



4>\ and 4>2 are complex scalar fields and A^ is an Abelian 
gauge field. The parameter r\ gives the symmetry break- 
ing scale. This Lagrangian possesses all the necessary 
ingredients to give rise to semilocal strings and so allows 



insight into more complicated models, such as D-term 
inflation in N = 1 or N = 2 supersymmetric (SUSY) 
models [HEH and D3/D7 brane inflation [1|. 

Since this model is just the Abelian Higgs model with 
an extra Higgs field, the semilocal strings arising from 
it can be seen as embedded Nielsen-Olesen strings [44j . 
The vacuum manifold is S , which is simply connected, 
and the existence and stability of the strings depend on 
energetic and dynamical causes, not on the topology of 
the vacuum. There is only one parameter governing the 
stability of the strings, namely, 1 = A/2e 2 : for j3 < 1 
they are stable (even quantum mechanically); for (3 > 1 
they are unstable. The BPS case (/3 = 1) is the transi- 
tional case, where the strings are neutrally stable [la. kg. 
There is a zero mode concerning the string width in the 
BPS case: all strings with widths that range from the 
width of an Abelian Higgs string to infinity are degen- 
erate in energy. But even though the string can have 
a very large (infinite) core, the magnetic flux is never- 
theless quantized. It is also worth noting that when two 
semilocal strings of any width collide, they intercommute 
[471 . l48l . 1491 ] and their widths may revert to the minimum 
width (i.e., the width of an Abelian Higgs string) 48]. As 
mentioned earlier, the BPS case arises naturally in, for 
example, D-term inflationary scenarios, In more generic 
cases, non-BPS cases will arise. We can anticipate some 
qualitative behaviour for the non-BPS case. For j3 > 1 
the semilocal string become unstable, and they will be- 
have like global textures. For j3 < 1, the lower the value 
of (3, the more semilocal string networks resemble Abelian 
Higgs string networks [4fJ, and their CMB predictions 
will tend to that of Abelian Higgs strings. 

One important difference between Abelian Higgs cos- 
mic strings and semilocal strings is that the latter are 
non-topological, so they can have ends. We therefore do 
not expect infinite semilocal strings to be formed by the 
Kibble mechanism in a cosmological phase transition. In- 
stead, a collection of segments forms. The evolution of 
the segments is very complicated: some will shrink to 
zero, others will join both ends to form loops, while some 
will join to nearby segments to make longer segments, 
and eventually form virtually infinite strings [4(| . At the 
BPS limit, the formation of longer segments is strongly 
suppressed, which led the authors of Refs. [19|,[2(| to con- 
jecture that no strings would be formed in such a model 
after inflation, and to claim that the D-term inflation 
without cosmic strings could take place. 

While the semilocal model can be reduced to the 
Abelian Higgs model by removing one complex scalar 
field, it can also be reduced to a global texture model by 
removing the gauge field (see Appendix) . A key question 
is how those removals modify the predicted anisotropy, 
keeping fundamental Lagrangian parameters fixed. An 
important derived parameter is the string mass per unit 
length in the Abelian Higgs and semilocal cases [4J, [50] : 



/! = 27T77 
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1 See for instance Ref. [41] for a comparison of the dynamics of 
the linear and non-linear a models. 



and the CMB anisotropies are proportional (Gfi) 2 , where 



G is the gravitational constant. We continue to use this 
measure even when talking about the texture model, but 
in that case it has no direct interpretation and is merely 
a measure of the energy scale of the model. It is known 
that textures give significantly lower anisotropies than 
strings for a given G/x, since textures decay much more 
quickly once inside the horizon. Our original precon- 
ception was that the anisotropies from semilocal strings 
would be of the same order as, but a little smaller than, 
those of textures. This expectation is because the gauge 
field cancels some of the field gradients present in the 
texture case and contributes little to the energy density 
[46j |. while in the BPS limit there are no long strings and 
hence no expectation of large anisotropies from string- 
like sources. Indeed, we will find that this expectation is 
confirmed by our detailed analysis. 



III. CMB CALCULATION METHOD 
A. Overview 

As discussed in the introduction, the method used to 
calculate CMB power spectra employed here is exactly 
that used for Abelian Higgs strings in Ref. [37| and we re- 
fer the interested reader to that article for greater depth. 
However we will briefly describe the fundamentals of the 
approach in this Section. 

Since observations reveal that the cosmological per- 
turbations are small (except on short scales in the cur- 
rent epoch) the defects can be taken to evolve in an 
expanding flat Friedmann-Lemaitre-Robertson-Walker 
universe. Inflation is taken to set up small primor- 
dial fluctuations, which then evolve passively under the 
Einstein-Boltzmann equations. However, the defects 
perturb the metric, via their energy-momentum tensor, 
and these perturbations also have to be taken into ac- 
count when computing CMB power spectra. But since 
the perturbations introduced by either inflation or defects 
are very small, any coupling between them is insignificant 
and the two perturbation sets may be calculated sepa- 
rately. The result is two statistically independent sets 
of CMB anisotropies, the power spectra of which simply 
add to give the total. 

Fundamentally the CMB power spectra are two-point 
correlation functions of the CMB temperature (or po- 
larization) fluctuations, and in the defect case they can 
be determined from the two-point unequal-time correla- 
tion functions (uetcs) of the defect energy-momentum 



tensor T_ [35j, IM |5J, \5c 
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Here k is the comoving wavevector, r and r' are two 
conformal times and ~ denotes a Fourier space quantity, 
as in the notation of Ref. [37| . 

However, the above 55 complex functions of 3 variables 
may be reduced to just 5 real functions of two variables, 



via the use of statistical isotropy, energy-momentum con- 
servation, and a property of defect networks known as 
scaling. The last property derives from the causal na- 
ture of the network: that field smoothing is limited by 
the horizon and, as such, defect networks often tend to 
an attractor regime in which statistical measures of their 
distribution are a function of a single quantity: the hori- 
zon size r. Assuming scaling and statistical isotropy, the 
uetcs can be written as [53, K33I 



UKXnu(k,T,r) = 
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4 CW„(fcT,fcr') (4) 



where V is the fiducial simulation volume and C is the 
scaling function of the uetc. Further, the functions C 
decay for large or small t/t' or for large ky/rr' and hence 
their measurement is only required over a fairly limited 
parameter space. 

Of the 5 independent scaling functions, 3 of 
them represent scalar degrees of freedom (Cf 1 (fcr, fcr'), 
Cf 2 (kr, fcr'),C'|2(fcT, kr')) while the remaining two repre- 
sent the vector and tensor degrees of freedom respectively 
(C v (fcr,fcr'), C T {kT,kT'j) [Hil. Unlike the inflation- 
ary case, the amplitude of the tensor perturbations is 
fixed and directly obtained from the simulations, and it 
is not possible to neglect the vector perturbations as they 
are continuously sourced by the presence of the defect 
network. 

The uetc scaling functions are then fed into a modi- 
fied version of the CMBeasy Boltzmann code [54| which 
computes the CMB temperature and polarization power 
spectra. However, in order for CMBeasy to use the 
uetc data, the functions C(kr, kr') must be first de- 
composed as: 



C(kT,kr') = ^A n c„(fcr) c n {kr') 



(5) 



which in the numerical case is tantamount to deter- 
mining the eigenvalues and eigenvectors of a real, sym- 
metric matrix. Note that Cf 2 is not symmetric, but 
Cfi (fcr, kr') = Cf 2 (kT , 1 kr) and so the M x M matrices 
representing the scalar degress of freedom (Cf 1; Cf 2 jCfi 
and Cf 2 ) are tiled together to yield a symmetric 2M x 2M 
matrix which is then decomposed this way. 

The modified CMBeasy then takes a single eigenvec- 
tor c„ as a source function and determines the corre- 
sponding contribution to the CMB power spectra, with 
summation over all such contributions yielding the to- 
tal. We found that for the semilocal case, the number of 
eigenvectors that should be included in the computation 
of the power spectra to achieve convergence was lower 
than that needed in the Abelian Higgs case; nevertheless, 
we use the same number of eigenvectors as in the Abelian 
Higgs case (see Table [TJ to minimize the differences in the 
computation. Also there are some choices about how to 
decompose the uetcs into the form of Eq. $5§ , and again 
we follow exactly the method of Ref. [37| . 



Note that scaling is broken at the radiation-matter 
transition and therefore the scaling UETC functions must 
be determined for the defects under both matter and 
radiation-domination. These are decomposed into eigen- 
vectors and fed into the modified CMBeasy code, which 
then interpolates between the radiation and matter re- 
sults to model the transition. For the transition into the 
A-dominated epoch, the sources simply decay and simu- 
lation is not required. In any case, perturbations sourced 
at such late times affect only the very largest scales. 



B. Semilocal simulations 

In order to calculate the scaling UETC functions for 
the semilocal string model above, we perform field 
theory simulations under both radiation- and matter- 
domination. However, since at times of importance for 
CMB calculations the string width is many orders of 
magnitude smaller than the horizon size, it is not pos- 
sible to perform simulations for such times while still 
resolving the string cores. Fortunately, scaling implies 
that this is not required and that the scaling uetc func- 
tions may be obtained from simulations of early times, 
with r ~ 10077" 1 , and then used to describe the late- 
time behaviour of the system. Note that although we 
perform simulations in a matter-dominated background 
for t ~ 100?y -1 , when the real universe was radiation- 
dominated, scaling implies that the network properties 
are the same at lOO^ -1 as they are in the true late-time 
matter era. 

However, performing the simulations required with 
current or near future computational technology is a con- 
siderable challenge. The (minimum) string width (r m i n ) 
is fixed in physical coordinates and, in the comoving co- 
ordinates required for the simulations, it varies as: 
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where a is the cosmic scale factor. It therefore shrinks 
as r _1 under a radiation- or as r~ 2 under a matter- 
dominated universe. If the strings are resolved at the 
end of the simulation, then they are enormous and larger 
than the horizon size for times only a little earlier, par- 
ticularly under matter domination. The system will not 
scale, or even contain significant regions where the fields 
lie on the vacuum manifold, if r is comparable to r m j n . 
Thus, the uetc data can only be collected over a limited 
region of the required parameter space. 

A simple means of circumventing the above problem 
is to promote the parameters A and e in the action to 
time-dependent variables [37| as: 



\ \ -2(l-s) 

e = e a- (1 - s) . 



Now the scale r m -, n varies as: 
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and variation of the action yields equations of motion as: 
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For the case of s = 1, the true dynamics are obtained, 
while for s — the factors of a responsible for the shrink- 
age of the (minimum) string width are removed from the 
equations of motion and r m i n becomes comoving. This is 
similar to methods of Refs. [55|, |5|| , where the undesir- 
able factors of a were directly removed from the dynami- 
cal equations, but here the parameter s can be varied and 
its effects measured. This is particularly important be- 
cause s 7^ 1 means that energy-momentum conservation 
is violated and it is T M „ that sources the perturbations. 
For the simulations under radiation-domination we can 
actually simulate the s = 1 case directly and compare 
in detail our intermediate results for s < 1 with the true 
case. However, one of the results of Ref. [37| was that the 
CMB results for Abelian Higgs strings are insensitive to 
s, and here we extend that result to semilocal strings. 

There is good evidence from previous simulations of 
both Abelian Higgs and semilocal strings that the net- 
work results are independent of the initial conditions once 
scaling is reached [33, [5a, [57], HI, [37j . The goal is there- 
fore to produce an initial configuration in a numerically 
feasible way that obeys the analogue of Gauss' law in 
the model and that will yield scaling fairly fast, in or- 
der to maximize the dynamical range of the simulation. 
Following again Ref. [37[ we set all temporal derivatives 
and gauge fields to zero, and the fields 4>\ and 02 are 
set in such a way that they lie in the vacuum manifold 
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with random phases in S 3 . The sys- 



tem is then evolved using the equations of motion ob- 
tained after discretizing the action via the standard Mo- 
riarty et al. [59j scheme and encoded in C++ using the 
LATfield Library [6y]. This scheme preserves the Gauss 
constraint during the evolution, even in the case of s < 1 , 
because our equations are derived from the discretized 
action (see Ref. 13711 1. 



IV. SIMULATION AND CMB POWER 
SPECTRUM RESULTS 

A. Numerical aspects and tests for scaling 

The numerical simulations were performed on the UK 
National Cosmology Supercomputer [6l|, where we sim- 
ulated five different realizations for both matter- and 
radiation-dominated epochs and for s — 0.0 and s = 0.3, 
as well as s = 1 for radiation alone. Of those, only the 
simulations with s — 0.3 were used to produce the CMB 




FIG. 1: Variation of r m i n and the parameters e and A for the 
s = 0.3 simulations, which mimic the one used in Ref. [371 ] . 
The subscript denotes the value at the end of the simulation. 



Parameter 



Value 



Lattice size N 

Lattice spacing Ax [rj^ 1 ] 

Time-step At [r/^ 1 ] 

Scalar coupling Ao 

Gauge coupling eo 

Initial conformal time n [v -1 ] 

Final conformal time r en( j [rj~ ] 

Initial s 

Final s 

Time of change in s value [rj ~ ] 



512 
0.5 
0.1 
2.0 
1.0 
1.0 
128 
-0.116 
0.3 
32 



T sim range of £ fit [r? -1 ] 
Dynamic range (r/r') ma x 
UETC sample matrix size 
UETC matrix sampling 
No. eigenvectors 



50 - 128 

~2 

512 

linear 

128 



TABLE I: Parameters used in production runs. The first 
part of the table lists the parameters of the simulations, with 
the second part listing parameters of the UETC method for 
calculating the CMB power spectrum. Note that the values 
are identical to those of Ref. [37]], except for the time range 
of the simulations. 



power spectra; the others were used to test the validity 
of our approach. We begin and end the simulation with 
the same value for r m j n ; that is, we initially set s to a 
negative value so that r m j n increases until a time r s , af- 
ter which s is positive so that r m j n shrinks during the 
primary part of the simulation (see Figfl]). This process 
accelerates the formation of vacuum regions. 

Table U shows the values of the parameters used for the 
production runs. Most of the parameters match those of 
Ref. [37|] ; furthermore, the coupling parameters were also 
varied as shown in Fig. [1] to replicate the behaviour of 
Ref. [33] . Because of the existence of zero modes in the 
semilocal system it is dangerous to run the simulation for 
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FIG. 2: Average value of £, as defined in Eq. (|10[) . for 
s = 0.0, 0.3, 1.0, in the radiation era. The average is over 
five different realizations for each value of s, and the shaded 
regions correspond to 1-er and 2-er deviations. The best-fit 
line for times r > 50 is also shown (dashed line). Note that 
the best-fit line is an excellent approximation, and that the 
differences between different values of s are minimal. 



longer than the half-crossing-time. In order to keep the 
dynamic range of Ref. (37| , here we had to begin taking 
UETC data at an earlier time (r = 50 rather than r = 64). 
One might wonder whether the system is already into the 
scaling regime by those early times, but we will show in 
this Section that in fact this is the case. 

One means to test for scaling is via the string length 
density, which should be proportional to r~ 2 , since the 
average length of string per horizon volume should be 
proportional to r and the horizon volume varies as r 3 . 
The detection of (topological) Abelian Higgs strings in 
a lattice field simulation is fairly straightforward, since 
one can track, for instance, zeros of the Higgs field or 
the winding of its phase around them. Alternatively one 
could use the average Lagrangian density (as in Ref. [371 ]) 
or the Too component of the energy-momentum tensor, 
which should also vary as t~ 2 . 

In fact, for (non-topological) semilocal strings those 
alternative methods must be employed: the zeros of the 
Higgs fields are not a good measure, since there is no need 
for the string core to be a zero of the Higgs field. Besides, 
semilocal strings in the Bogomol'nyi limit ((3 = 1) are ex- 
pected to form short segments that either collapse very 
quickly or become fatter, diluting away the core of the 
string. That is why we used a different measure to quan- 
tify the scaling of the semilocal BPS string system. 

In this work, we use the Too component of the energy- 
momentum tensor to define our length measure £ as: 



i 
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(10) 



and use it to check the scaling of the system. Figure [2] 




FIG. 3: The dependence of the 5 ETCs upon s in the radiation era at r s i m = 1287? -1 . The shaded regions show the 1-cr and 
2-cr uncertainties for the s = 0.3 case, while the lines indicate the s = 1.0 (solid), s = 0.0 (dot-dash) results, which clearly lie 
within the statistical uncertainties. 



shows the result of the simulations for the radiation era 
for s — 0.0, 0.3, 1.0, as well as the best fit for times 
t > 50. First of all, note that the differences with respect 
to varying the value of s are minimal. It is also clear from 
these figures that after the initial non-scaling period, £ 
follows a linear behaviour 



£ oc (t-T£ =0 ). 



(11) 



There is nothing fundamental about the offset r^ = o, as 
argued for Abelian Higgs case in Ref. j37|; it is a con- 
sequence of the non-scaling initial period of the simula- 
tions and is negligible in the late-time limit of interest 
for CMB calculations. What is fundamental is the value 
of the slope of £, which shows little dependence upon 
s: 0.10 ± 0.01 for s = 0.0, 0.10 ± 0.02 for s = 0.3 and 
0.09 ± 0.02 for s = 1.0. The difference between the three 
plots in Figure [2] is due almost entirely to differing values 

of r ?=0 - 

However, given that the system scales not with r but 
with t— T^ = Q, then following Ref. [37] we perform a rescal- 
ing of the measured UETC functions C(kr, kr') in accor- 
dance with: 



T - T ?=0 - 



(12) 



These functions then provide a more important test of 
scaling for the present work: after the rescaling, do they 
show negligible absolute time dependence? Figures [3] and 
2] show the equal-time correlator (etc) scaling functions, 
that is when t = t' . These are the most important part 




FIG. 4: The equal-time scaling function Cfi(kr, kr) averaged 
over five realizations for s — 0.3 in the radiation era. The dif- 
ferent lines correspond to roughly uniformly spaced r values 
in the range 50r)~ 1 < r < 12877 _1 . The shaded regions show 
the 1-cr and 2-cr uncertainties in the mean indicated for the 
latest time (solid). For visualization purposes, the mean offset 
across the five realizations is used, whereas the actual CMB 
calculations use independent offsets for each realization. 



of the UETCs, and the most suitable for visualization. 
The different lines in Fig. [4] correspond to the etc func- 
tion averaged over five realizations, for different times; 
the shaded regions correspond to the 1-cr and 2-cr un- 
certainties in the mean at the last time step (when the 
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FIG. 5: CMB temperature power spectrum from semilocal 
string simulations. The solid black line is the total prediction, 
which is obtained by adding the contribution of scalar (dash- 
dot red), vector (thin blue) and tensor (dashed green) modes. 



system is closest to scaling). Clearly, the different lines 
lie within the statistical uncertainties of the simulation, 
and there is no absolute dependence of the etc on time. 
As the uetcs are the only quantities entering the CMB 
power spectrum calculations, we are confident that the 
earlier sampling does not affect the results. 

We can also make use of the etcs to show that the de- 
pendence of the simulations with s is quite mild. Figure [3] 
shows the etcs for the last time step in the simulation 
for s = 0.0, 0.3, 1.0 for the five different uetcs need for 
the CMB calculation. In all cases the dependence on s 
lies well within the statistical uncertainties. 
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FIG. 6: CMB temperature power spectrum from Abelian 
Higgs strings (left) and texture (right) simulations. The solid 
black line is the total prediction, and its constituents: scalar 
(dash-dot red), vector (thin blue) and tensor (dashed green). 
Note the difference in the scale between the two graphs. 



Temperature power spectra from semilocal 
strings 



After checking that our numerical approach for simu- 
lating semilocal strings is satisfactory by monitoring the 
influence of the parameter s in the simulations and con- 
firming the early onset of scaling, the resulting uetcs 
were fed into the modified CMBeasy Boltzmann code 
as explained in Section IlIII The perturbations obtained 
from the strings and from inflation are calculated sep- 
arately and added in quadrature - the tiny interaction 
between the two contributions can be safely neglected. 
Besides, it is numerically quite slow to recompute the 
string spectra for many different values of the cosmologi- 
cal parameters. For these reasons, the string spectra were 
computed for a fixed value of the cosmological parameters 



(the central values given by cosmological experiments"). 
We then tested at the end that the string spectra do not 
change significantly when the cosmological parameters 
are varied within the limits imposed by the CMB data. 
Note that this is the only point where the cosmological 
parameters are involved in the calculation of the CMB 
spectra. 

Figure [5] shows the CMB temperature power spectrum 
from semilocal cosmic strings, divided into scalar, vector 
and tensor components. For comparison, we plot the cor- 
responding power spectra from Abelian Higgs strings and 
textures in Fig. O The power spectra for Abelian Higgs 



The calculation of the string spectra was done for h = 0.72 |62j . 
Q h h 2 = 0.0214 [H, H A = 0.75 [U. We also assumed spatial 
flatness and the optical depth to the last-scattering surface r = 
0.1. 
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FIG. 7: The temperature power spectra from semilocal strings 
(thick) and textures (thin), showing the scalar (dash-dot red), 
vector (dashed blue) and tensor (lower solid green) compo- 
nents. 
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6000 



Normalization of different defects to match the 
= 10 multipole value. 
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strings are taken from Ref. [33], 3 whereas the texture 
spectra have been calculated by using new simulations of 
the texture linear-cr model obtained by setting to zero the 
gauge fields in the semilocal model Eq. fl} (see Appendix 
for details on the texture calculations). 

The first thing to notice from these graphs is that, 
at fixed G/i, the Abelian Higgs strings give significantly 
higher anisotropics than the other two by a factor of 
about 5. Semilocal strings and textures are quite sim- 
ilar, with the former being slightly lower. This is in good 
accord with the intuition that semilocal strings ought to 
be similar to textures, but with some of the anisotropics 
cancelled by the gauge field. Ultimately, this difference 
in anisotropy amplitude will result in a weakening of ob- 
servational constraints on G/i for semilocal and texture 
as compared to Abelian Higgs strings. 

Going beyond the amplitude, it is clear that the semilo- 
cal temperature power spectra shape shares properties 
with both the Abelian Higgs and the texture prediction. 
It is close to the texture prediction in, e.g., the level of 
the power of the anisotropies at fixed G/i, although it 
does not have the small oscillations. On the other hand, 
it peaks for larger values of £ than textures, more like 
Abelian Higgs strings. 

Figure[7]shows a direct comparison between the semilo- 



3 Refs. Il5l |37H employed a code in which a bug has been discov- 
ered, and this had a small effect in Ref. [LUl since it used their 
results directly (see the respective errata) . Here we have used 
the corrected power spectra from Refs. [l5l |37I | and quote the 
corrected results from Ref. |13|| . 



FIG. 8: The temperature power spectra for Abelian Higgs 
strings, semilocal strings, textures and the best fit inflation- 
ary model. The poi nts represent the three-year WMAP ex- 
perimental data [65}. The curves for defects are normalized 
to match the data at £ = 10. 



cal strings prediction and the texture prediction for the 
temperature power spectrum. As discussed before, the 
texture model can be obtained by setting to zero the 
gauge coupling constant in the semilocal model, so the 
differences appearing in the curves arise from coupling 
the scalar fields to the gauge fields. The overall shape of 
the curves for each model is roughly the same, with that 
for textures being somewhat higher and peaking earlier 
than that for the semilocal strings. If the Abelian Higgs 
model prediction was plotted in the same graph, it would 
lie way above this curve, since the overall scale of the 
spectrum is much higher. 

However, bear in mind that the normalization of the 
spectra is a free parameter, and one would like to com- 
pare the relative shapes of the predictions to understand 
their detailed observational implications. For example, it 
is customary to encode the information about the overall 
power coming from the defects by G/xin. The normal- 
ization of the strings is proportional to (G^i) 2 . Thus, 
Gfiio is the value of the normalization by which the de- 
fect spectra matches the WMAP data [65| at multipole 
£ — 10, shown in Table HT1 for each type of defect. This 
is merely a tool to compare different predictions. Fig- 
ure [8] shows pictorially the different defect predictions 
with these normalizations, together with the experimen- 
tal data and the inflationary best-fit model. Note that 
for textures it is not natural to talk about "energy per 



unit length" , but we adopt the same normalization strat- 
egy Eq. @ by analogy with the strings, in order to ease 
the comparison between models. Different normalization 
schemes are discussed in the Appendix. 



C. Polarization power spectra from semilocal 
strings 

The polarization power spectra are shown in Fig. [5J de- 
composed as usual into the ee , te and bb components 
|6a |. The curves are normalized so they have equal tt 
power at £ — 10 (and match the WMAP measurement 
on this scale). 4 Note that, as in Ref. [l5| . we do not in- 
clude the coupling between inflation and defect perturba- 
tions as a result of gravitational lensing at late times. We 
do not, however, expect this effect to be significant: the 
gravitational lensing of the inflationary scalar ee spec- 
trum by the defect perturbations is likely to be insignifi- 
cant relative to that produced by the inflationary pertur- 
bations, since defects are sub-dominant. It is also likely 
to be negligible relative to the B-mode from defects since 
defects still produce a strong B-mode signal even when 
their sub-dominance is taken into account. The lensing of 
the defect ee and BB spectra by the inflationary pertur- 
bations should be negligible also since defects contribute 
equally strongly to these two spectra. 

With the temperature spectra all normalized to the 
same value at £ — 10, the polarization spectra are all 
of similar magnitude and shape (if normalized to the 
same G/i, the Abelian Higgs contribution would be much 
higher). In the ee case, all three defect models peak 
at the same value of £ ~ 1000, with the semilocal com- 
ponent having slightly less pronounced oscillations, and 
the secondary peak at low £ coming from reionization 
follows the following trend: it peaks at the lowest value 
of £ for textures, next for semilocal strings, and highest 
for Abelian Higgs strings. The same trend can be seen 
for the \ow-£ features of the te spectra. The high-£ be- 
haviour of the te is roughly the same for all three cases, 
showing large oscillations that eventually die out. Note, 
however, that the oscillations for the semilocal case are 
the smallest, and there is no alternation between anticor- 
relation and correlation due to the combined effect of the 
scalar-vector-tensor contributions which have different 
phases. 

The BB power spectrum is arguably the most inter- 
esting prediction from defects [lj, [l5|, [l8|, [3511: inflation 
contributes weakly to the bb polarization [67] . and the 
signal predicted from defects might be detected in the 
near future, providing a smoking gun for the existence 
of defects in cosmology. In Fig. [5] it can be seen that 






O 




FIG. 9: Polarization power spectra for semilocal strings (solid 
red), compared to Abelian Higgs strings (dashed black) and 
textures (dot-dash blue). The top figure is the te power 
spectra, the ee is in the middle, and BB at the bottom. All 
the curves are normalized by making the temperature spec- 
tra match the WMAP £ = 10 value, i.e., using the values of 
Gfi — G/iio as in table [HI In the te figure, the thick lines 
correspond to correlations and the thin lines to anticorrela- 
tions. The results plotted here correspond to a value t — 0.1 
(see footnote [5]). 



the shape of the BB power spectra is similar for semilocal 
strings, textures and Abelian Higgs string. Both peaks 
follow the aforementioned trend, with textures peaking 
for lowest £ ~ 240, then semilocal strings £ ~ 360, and fi- 
nally Abelian Higgs strings £ ~ 590; furthermore, the BB 
power spectrum for textures decays first, then that for 
semilocal strings, and the Abelian Higgs string bb power 
spectrum decays at the highest £. 



CONSTRAINTS FROM CURRENT CMB 
DATA 



4 See Fig. 1121 in Section [V] for a figure with the normalization 
of each case set to match the respective 95% upper limits from 
current data. 



Clearly, a string component at the G/xio level shown 
in Table [H] is ruled out by the experiments. There- 
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Parameter 


PL 


PL±SL 


PL±AH 


PL±TX 


/io 





0.14 ±0.06 


0.09 ±0.05 


0.18 ±0.08 


n s 


0.96 ±0.02 


1.02 ± 0.04 


1.00 ± 0.03 


1.04 ±0.05 


h 


0.74 ±0.03 


0.85 ± 0.06 


0.84 ± 0.06 


0.87 ±0.08 


n b h 2 


0.0223 ± 0.0007 


0.0263 ± 0.0021 


0.0258 ± 0.0020 


0.0264 ± 0.0026 


n m h 2 


0.126 ±0.007 


0.122 ±0.007 


0.122 ±0.007 


0.118 ±0.008 


ln(10 10 A 2 s ) 


3.09 ± 0.06 


2.96 ± 0.08 


2.99 ± 0.08 


2.94 ±0.09 


r 


0.09 ± 0.03 


0.11 ±0.04 


0.11 ±0.04 


0.12 ±0.04 



TABLE III: Mean and standard deviations of the power-law model parameters, using CMB data only. 
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FIG. 10: The 2-dimensional constraints from the MCMC 
chains using CMB only data (bottom), and using 
CMB±BBN±HKP data (top). 



fore, in order to quantify the allowed level for a string 
component, we introduce a parameter /io [131 ] : the frac- 
tional contribution of strings to the temperature power 
spectra at multipole £ — 10. We then use a modified 
version of CosmoMC [68| to perform a multi-parameter 
Markov Chain Monte Carlo (MCMC) likelihood analysis 
for CMB data (WMAR ACBARJ300MERANG, CBI 
and VSA projects [H, M, M, M, IB, Izl ) when semilocal 
strings or textures are included. We also reproduce the 
results for Abelian Higgs strings from [13j for compari- 
son. The likelihood analysis is performed for the usual 
six-parameter Power-Law (PL) model (Vt\,h 2 , h, n s , r, 
O m /i 2 , A s ), plus the new extra /io- We vary only the 
normalization of the string (or texture) power spectra 
and keep the form of them fixed, following Refs. |3.[l3|. 
This approach is justified since the strings (or textures) 
contribute only a small fraction to the total power, and 
the changes in the cosmological parameters are not large 
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FIG. 11: The 1-dimensional marginalized likelihood plots for 
the fraction of semilocal strings /io allowed at multipole t = 
10 from the MCMC chains using CMB only data (dark), and 
using CMB±BBN±HKP data (light). 



enough that the changes in the form of the string or tex- 
ture spectra are relevant. The best- fit values of the pa- 
rameters can be found in Table IIIII for PL alone, and 
for PL with semilocal (SL) strings, Abelian Higgs (AH) 
strings and textures (TX). 

Table Hill shows that, as in previous results [9|, uM, the 
degeneracies of the cosmological parameters with respect 
to the CMB data allow for a rather high value of fio — 
0.14±0.06, which corresponds to Gfi = [1.9±0.4] x IO" 6 . 
In order to accommodate the defect contribution, other 
parameters (most notably £lbh 2 , h, and n s ) get shifted 
to higher values, as was the case for the Abelian Higgs 
strings 13]. 

Figure [JO] shows the marginalized 2D likelihood distri- 
butions for the parameters £lt,h 2 , h, and n s versus /io- 
It is clear in the figure that there is ample space for a 
noticable value of fio. In fact, it can be seen in the 1- 
D marginalized likelihood plot (Fig. [TT]) that the CMB 
data alone prefers to have semilocal strings at a 2-cr level. 

However, as mentioned before (see Table IHT[) . the best- 
fit values of U^,h 2 and h are high compared to the con- 
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cordance model: Non-CMB experiments for those two 
quantities yield h — 0.72 ± 0.8 (the Hubble Key Project 
HKP 62]) and tt b h 2 = 0.0214 ± 0.0020 (measurements 
of deuterium abundance in high-rcdshift gas clouds to- 
gether with big bang nucleosynthesis BBN |63|). If we 
include these two non-CMB experimental values into our 
MCMC scheme (as gaussian likelihoods), we confirm 
that the new set of data is more constraining also for 
semilocal strings (as seen in Figs. [TOl and ITTI) . However, 
the 95% confidence level upper bound on f w for semilocal 
strings remains rather high, going from Gp, < 2.6 x 10" 6 
(/io < 0.25) for the CMB only case to Gp < 2.0 x 10" 6 
(/io < 0.14). 

We may now use these upper bounds to compare the 
strengths of the power spectra of different models. Ta- 
ble |IV] summarises the values of /io and G\x we obtain for 
our analysis of semilocal strings and textures, together 
with that for Abelian Higgs strings from Ref. 13]. The 
parameter /io depends mainly on the form of the de- 
fect power spectra, and we see that the data allow for a 
higher contribution from semilocal strings (and textures) 
than from Abelian Higgs strings. Besides, the values 
of G/iio are also higher for semilocal strings (and tex- 
tures) than for the Abelian Higgs strings. Therefore, the 
expectation that the constraints on the Abelian Higgs 
string tensions would be alleviated by semilocal strings 
is confirmed: the 95% upper bound limit for semilo- 
cal strings is Gp < 2.0 x IO" 6 for CMB+BBN+HKP 
(Gfi < 2.6 x 10" 6 for CMB only) whereas for Abelian 
Higgs string it is Gp < 0.9 x IO" 6 for CMB+BBN+HKP 
(G/i < 1.1 x 10" 6 for CMB only). 

Figure [T2] shows the temperature power spectrum 
and the BB polarization spectrum for semilocal strings, 
Abelian Higgs strings and textures, normalized by using 
the upper bounds obtained from performing the MCMC 
likelihood analysis for both CMB and non-CMB (BBN 
and HKP) data. In the temperature power spectrum fig- 
ure it can be seen that due to the different shapes of the 
curves, even though the textures are allowed the high- 
est level of power at £ = 10, the semilocal and Abelian 
Higgs strings go up to higher peaks. The signal is high- 
est for textures for £ < 60, for semilocal strings for 
60 < £ < 250 and then for Abelian Higgs for 250 < L On 
the other hand, for the bb polarization, the normaliza- 
tion makes the trend of the power coming from the three 
cases just the reverse of the temperature case: textures 
have the highest peak, followed by semilocal strings and 
then Abelian Higgs strings. Needless to say, the normal- 
ization does not change the value of I at which each case 
peaks, and we find that textures peak first, then semilo- 
cal strings, and then Abelian Higgs strings. 
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FIG. 12: Temperature power spectra (top) and bb polariza- 
tion power spectra (bottom) for semilocal strings (solid red), 
compared to Abelian Higgs strings (dashed black) and tex- 
tures (dot-dash blue). The curves are normalized using the 
95% confidence level upper bound from a MCMC analysis 
using CMB , BBN and HKP data. 



Model 


fw 


CMB only 
Gfj, 


n s 


CMB+BBN+HKP 
/io Gfj, n B 


SL 
TX 
AH 


0.25 
0.33 
0.17 


2.6 x 10"" 
2.5 x 10"" 
1.1 x 10"" 


1.09 
1.14 
1.06 


0.14 2.0 x 10~" 1.01 
0.16 1.8 x 10"" 1.02 
0.10 0.9 x 10"" 1.00 



TABLE IV: 95% upper bounds on f w , G\x and n 3 for PL+X, 
using CMB only, and using CMB+BBN+HKP. 



out 5 (see table HV)) . Therefore, when comparing differ- 
ent cosmological models, one should also consider the 
(scale-invariant) Harrison-Zel'dovich (HZ) model with 5 
parameters (fib^ 2 , h, r, f2 m /i 2 , A,) and a fixed n s = 1. 

Table [V] shows the goodness of fit for the concordance 
PL model, as well as for the HZ model, when Abelian 
Higgs strings, semilocal strings and textures are added; 
for both CMB only and CMB+BBN+HKP data. In 
most cases, the inclusion of the defects makes the fit bet- 
ter, with semilocal strings obtaining marginally the best 
fit. This is not too surprising for the cases where defects 



It is worth mentioning that, in all cases studied, the 
value for the spectral index n s is higher than the one ob- 
tained for cases without defects; and n s — 1 is not ruled 



See Ref. [7J| for a discussion about n B in light of WMAP3 from 
a Bayesian model selection point of view. 
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CMB only 


CMB±BBN±HKP 


Model 


Ax 2 


ln(E) 


Ax 2 


HE) 


HZ 


6.6 


-0.6 ±0.1 


8.4 


-1.4 ±0.1 


PL 














HZ+SL 


-4.4 


2.5 ±0.2 


-0.7 


1.0 ±0.3 


PL+SL 


-4.5 


1.0 ±0.3 


-1.9 


-1.0 ±0.1 


HZ+AH 


-3.5 


1.8 ±0.2 


1.4 


-0.7 ±0.2 


PL+AH 


-3.5 


-0.1 ±0.2 


-1.4 


-1.8 ±0.2 


HZ+TX 


-3.2 


2.0 ±0.3 


-0.5 


0.3 ±0.3 


PL+TX 


-3.4 


1.2 ±0.3 


-1.2 


-1.2 ±0.1 



TABLE V: Goodness of fit and Bayesian Evidence relative to 
the concordance model PL. 



are added to the PL model, because we are adding an ex- 
tra parameter to the PL model, and of course obtaining 
a better fit. But the models with HZ ± defects have the 
same number of parameters as the PL model (6 parame- 
ters), and surprisingly, the fit is better when considering 
CMB data only. 

In order to compare different models, we use Bayesian 



evidence values [75[. The evidences are calculated us- 
ing the Savage-Dickey method [76|, [77j with flat priors 
of < /io < 1 and 0.75 < n s < 1.25 (Table EJ. In 
the case of CMB data only, adding a defect contribution 
to the PL model makes little change to the evidence, 
whereas HZ±defects models have a significant higher ev- 
idence (and semilocal strings are marginally better than 
textures or Abelian Higgs strings). Therefore, taken at 
face value, HZ±SL strings is the six-parameter model 
which fits the data best, though the others are not con- 
vincingly worse. With the non-CMB data added, the 
evidences get reduced, but not to the level to conclude 
that a string contribution is ruled out. 

This type of comparison will have some dependence 
on the choice of prior used for the defect component (the 
other priors are less important as they are shared by the 
models). Our imposition of a prior on /io is just one 
out of several plausible options. A prior on fi might be 
more appropriate but is hard to formulate; for instance 
a logarithmic (Jeffreys) prior on /i would put most of the 
prior weight on models with negligible defect contribu- 
tion and hence leave the data unable to distinguish them 
(the same situation occurs with the inflationary tensor 
component - see Ref. [7J|)- ^ n an y event, we can safely 
say that present data are not good enough to exclude a 
defect contribution, and leave open the tantalising possi- 
bility that such a contribution may be a significant one. 



VI. DISCUSSION 

In this article we report on the first ever semilocal 
string simulations to predict CMB power spectra. This 
has been possible because of recent improvements in sim- 
ulating field theoretical defects, as there is no Nambu- 
Goto type approximation for semilocal strings. The case 



of Abelian Higgs model strings is much better known, 
so we compared the semilocal string predictions to the 
Abelian Higgs string ones, trying to minimize the dif- 
ferences in the approaches. Since semilocal strings were 
expected to share properties with both an Abelian Higgs 
string prediction and a texture prediction, we also pro- 
duced a set of simulations for texture- type defects to com- 
pare all three defect contributions. 

Our first conclusion is that the amplitude of the spectra 
for the semilocal strings is indeed lower than in Abelian 
Higgs models (for a given energy scale); therefore, the 
normalization at I = 10 to match the WMAP value gives 
a higher value for semilocal strings G/iio = 5.3 x 10~ 6 
than for Abelian Higgs strings Gfiio = 2.7 x 10~ 6 . More- 
over, the form of the power spectra for semilocal strings 
allows for their fractional contribution to the tempera- 
ture power spectra at multipole £ = 10 to be higher than 
the one for Abelian Higgs strings: the 95% confidence 
upper bounds are /io < 0.25 (/io < 0.14) for semilo- 
cal strings for CMB only (CMB±BBN±HKP); whereas 
/io < 0.17 (/io < 0.10) for Abelian Higgs strings. 

These two ingredients contribute to allow a higher 
value of G/j, for semilocal strings when CMB data (and 
also when HKP and BBN data) are taking into ac- 
count: Gfi < 2.6 x IO" 6 {Gfi < 2.0 x IO" 6 ) for semilo- 
cal strings for CMB only (CMB±BBN±HKP); whereas 
Gfi < 1.1 x IO" 6 (Gfi < 0.9 x IO" 6 ) for Abelian Higgs 
strings. Thus, if a promising high-energy physics infla- 
tion model faces a problem of creating Abelian Higgs 
cosmic strings at a scale that it is too high for the CMB, 
a possible resolution could be to extend the model to 
predict semilocal strings instead. This weakens the Gfi 
constraint by a factor of nearly three, which may give the 
theory extra breathing space. 

As in the case of textures [9| and local strings [ll|, 
[l3 |. the best fit for the cosmological parameters when 
taking into account semilocal strings gives a rather high 
spectral index n s , comparing to the concordance model. 
For CMB only data, the preferred value of n s is higher 
than 1 (n a = 1.02±0.04), whereas if BBN and HKP data 
are included n s = 0.98 ± 0.02. This is in contrast to the 
Abelian Higgs data, where the CMB data alone predict 
n s = 1; but the result is very similar when BBN and 
HKP data are included [13]. This result also alleviates 
pressure on high-energy physics models (such as D-term 
and F-term hybrid inflation) which predict defects and 
an n s close to unity. 

As n s = 1 is close to our best-fit values, we studied 
how well a combination of Harrison-ZePdovich plus de- 
fects would fit the data, comparing it to the concordance 
model. Using goodness of fit and Bayesian Evidence cri- 
teria, we find that models with defects are not ruled 
out (even when BBN and HKP data are included), and 
should be considered as competitive models. In fact, tak- 
ing into account only CMB data, there is significant ev- 
idence in favour of a Harrison-Zel'dovich plus semilocal 
strings model. 

We have not included large-scale structure datasets 
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into our analysis for two reasons. Firstly, perturbations 
from defects are non-linear and non-Gaussian as soon 
as they are created Q, and so it is not clear that the 
commonly- used model [78] for non-linear corrections to 
the galaxy power spectrum applies. Secondly, recent 
analyses [73, [SO] show scale-dependent bias even on rela- 
tively large scales (0.05 < /c/ZiMpc" 1 < 0.15), and make 
clear that more work is needed before the galaxy power 
spectrum can be used for robust estimates of cosmologi- 
cal parameters. 

The prediction for polarization power spectra from 
semilocal strings is similar to that of Abelian Higgs 
strings. As a matter of fact, the height of the bb polar- 
ization peak (normalized using the 95% confidence level 
upper bound on /io) is of the same order of magnitude 
for Abelian Higgs strings, semilocal strings and textures, 
though the position at which the maximum occurs vary 
from case to case. A detection of a bb signal on these 
scales could be an indicator of the existence of cosmic 
defects [i|. 

In the present work we have studied only BPS semilo- 
cal strings, but we will in a future work extend the anal- 
ysis to lower values of (3, which would arise in F-term su- 
persymmetric models (or even in P-term models [8UI82J). 
The behaviour of semilocal strings does depend strongly 
of (3 [40] : the lower the value of (3 the more the semilocal 
strings behave as Abelian Higgs strings. Tracking how 
the CMB predictions change with respect to (3 will give 
us insight into how the different components in a string 
network contribute to the creation of anisotropies. 

There is also a whole new class of strings coming from 
string inflation: (p, g)-strings |g, 0, HI- Even though 
there is some work about the networks of these defects 
[83L \84 . l85l l8q ] , basic properties, such as scaling, are still 
not really understood. Provided some sound CMB pre- 
dictions from these strings are calculated, it would be 
interesting to compare the results with those of Abelian 
Higgs or semilocal strings. 

Ultimately, one would like to be able to discriminate 
between all the different defect types, be it Abelian Higgs, 
semilocal, textures of (p, g)-strings. As we have shown 
here, the predictions do look similar, but there are differ- 
ences, e.g., the aforementioned position of the BB peak. 
The question is whether these differences are strong 
enough, or whether future data will be good enough, to 
pinpoint whether defects do appear and if so, which type. 
In a future publication we will address this issue. 



00042) and FPA2005-04823 (J.U.). The simulations 
were performed on COSMOS, the UK National Cosmol- 
ogy Supercomputer, supported by SGI, Intel, HEFCE 
and PPARC. 



APPENDIX A: TEXTURE SIMULATIONS 

In this work we have compared the CMB predictions 
to those from the Abelian Higgs model and texture mod- 
els. For the latter we have performed a new set of 
simulations for the linear-cr global texture model com- 
paring the results additionally with previous ones from 

Refs. [1, m, m, H3] . 

We exploited the fact that a texture model can be ob- 
tained directly from Lagrangian Eq. ([1]) by setting the 
gauge field to zero, obtaining 



C = 



A 



l^il' + l^^l'-^^il' + l^l 2 -^ 2 )' 



E(W 2 -z E^ 



2 2 

V 



(Al) 



where ipi, (i = 1...4) are real scalar fields and <f>\ = ipi + 

ilj} 2 , <f>2 = 4>3 +ilp4- 




multipole moment: I 



FIG. 13: The temperature power spectra from linear a tex- 
tures (thick) obtained and used in this article; together with 
the spectra from a non- linear a model (thin) from [371 ]. The 
different components are scalar (dash-dot red), vector (dashed 
blue) and tensor (lower solid green). 
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This is the Lagrangian we have simulated to obtain the 
figures used in the main body if this article. Note that 
it is not canonically normalized with respect to the real 
scalar fields, which becomes clear when compared to the 
non-linear a model: 






$ 



(A2) 
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where in this case A is a Lagrange multiplier fixing the 
fields onto the vacuum manifold. 

The normalization of the textures is usually given for 
the Lagrangian canonically normalized for scalar fields as 



e = 4ttGk 2 . 



(A3) 



Comparing both Lagrangians (|Alj) and (|A2j) . it can be 
seen that k 2 = 2r/ 2 , so the normalization scheme for tex- 
tures used throughout this work Eq. ([2]) relates to that 
of Ref . [13 by 



4Gft 



(A4) 



Note that Ref. [35( uses another definition of e, bigger by 
a factor 2ir. 



The procedure followed to obtain the CMB predictions 
from textures was exactly the same as for the semilocal 
case (and for the Abelian Higgs case) . The results can be 
seen in Fig. [13 Earlier work on global textures [4X|. l88i l89jj 
represented important steps towards the formalism we 
are using here. However, they neglected decoherence and 
so overemphasized the presence of the acoustic peaks. 

We also compared the texture predictions we obtain 
in our simulations with those of a non-linear a model by 
some of the authors of this article [33] in Fig. Q2J Those 
simulations were performed using different parameters; 
e.g., smaller lattices (256 3 ) and smaller number of eigen- 
vectors used to obtain the CMB power spectra. In any 
case, the agreement with our linear-cr model is very good, 
which makes our comparison of texture with semilocal 
and Abelian Higgs strings more widely valid. 
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